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SUMMARY 

On  9  June  1966  three  satellites  were  launched  into  polar  orbits  with 
perigee  heights  near  180  km  and  apogees  near  36OO  km.  They  were  Secor  6 
(l 966-51 B),  a  rectangular  box  of  mass  17  kg  and  length  0*33  m;  ORS  2 

(l 966-51 C),  a  regular  octahedron  of  side  0*23  m  and  mass  5  kgj  and  the  final- 

stage  Agena  rock*  t  (I966-5IA),  of  length  at  least  6  m.  Analysis  of  the  exten¬ 
sive  orbital  data  now  available  on  these  three  satellites  gives  260  values  of 
air  density  at  heights  near  1 90  km.  Comparison  of  these  values  provides  a 

severe  test  of  the  accuracy  of  the  orbital  data  and  methods  of  analysis,  and 

shows  that  the  accuracy  is  at  least  as  good  as  has  previously  been  estimated. 
The  comparison  also  reveals  that  the  effective  cross-sectional,  area  of  Secor  6 
oscillated  with  an  amplitude  of  about  ^Ofo,  with  minima  when  perigee  was  near 
the  equator  and  maxima  at  latitudes  near  70°.  As  a  consequence,  the  amplitude 
of  the  semi-annual  variation  in  air  density  previously  obtained  from  Secor  6 
should  be  reduced  from  1 *45  to  1 *40.  Other  findings  are  that  the  effective 
cross-sectional  areas  of  the  satellites  never  approached  their  extreme 
possible  values,  and  that  the  effective  area  of  I966-5IA  changed  by  less  than 
\C$  from  week  to  week. 
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INTRODUCTION 


The  three  satellites  1 966-51  A,  B  and  C,  launched  on  9  June  1 966,  have 

yielded  valuable  information  about  the  upper  atmosphere.  All  three  objects 

were  in  similar  polar  orbits,  with  perigees  near  180  km  and  apogees  initially 

near  3600  km.  The  decrease  in  inclination  of  all  three  orbits  was  used  to 

-] 

determine  the  rotational  speed  of  the  atmosphere  ;  the  orbit  of  I966-5IB 
2  3 

(Secor  6)  was  used  *  in  a  detailed  study  of  atmospheric  density;  and  the 
orbit  of  1966-51C  (ORS  2)  was  calculated  at  R.A.E.  over  a  short  time  interval, 
and  the  variations  in  density  analysed^. 

The  orbit  of  1 966-51 A  (Secor  6  rocket)  has  recently  been  determined^  at 
R.  A.  E. ,  mainly  from  visual  observations.  In  the  present  paper,  the  values  of 
orbital  acceleration  given  in  Ref. 5,  together  with  those  derived  from  Space- 
track  bulletins  for  1966-51A  and  C,  have  been  used  to  obtain  values  of  air 
density  for  comparison  with  the  results  of  Ref. 2. 

The  comparison  has  three  main  purposes.  First,  it  provides  a  severe 
test  of  the  accuracy  and  reliability  of  determining  density  from  Spacetrack 
orbital  data:  if  the  methods  or  the  data  are  unreliable,  the  densities 

2 

obtained  from  1 966-51 B  and  C  will  differ  significantly.  Second,  the  results 
from  1 966-51 B  indicated  an  important  semi-annual  variation  in  air  density,  by 
a  factor  of  1 *45  at  1 91  km:  this  is  the  only  such  evidence  available  for 
1966-7  at  heights  below  200  km,  and  the  value  needs  to  be  improved,  if  possible, 
because  it  has  a  crucial  effect  on  the  choice  of  upper-atmosphere  models^. 

Third,  the  area/mass  ratios  for  1 966-51 B  and  C  are  known,  with  errors  of  about 
±1C$,  and  the  results  can  be  compared  with  those  from  I966-5IA,  for  which  the 
size,  shape  and  mass  are  unknown,  to  find  the  area/mass  ratio  of  51A  and 
its  variations,  if  any. 

2  CALCULATION  OF  ATMOSPHERIC  DENSITY 

2.1  Orbital  decay  rate 

The  method  used  to  obtain  values  of  density  is  generally  the  same  as 
that  described  in  Ref. 2;  i.e.  the  Spacetrack  values  of  mean  motion  n  were 
differenced  to  give  values  An,  which  were  then  divided  by  the  time  interval 
At  to  obtain  values  of  h  =  An/At,  which  were  assumed  to  apply  at  the  mid¬ 
point  of  the  time  interval.  Figs. 1-3  show  the  values  of  h  for  1966-51B,  C 
and  A  respectively,  from  Spacetrack  and  other  source's. 
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The  only  difference  in  treatment  here  is  that  up  to  MJD  39347*8 
(1966  Aug  10)  the  values  of  n  given  by  the  R.A.E.  orbit  of  1 966-5'*  A  were 
used  .instead  of  Spacetrack  values  of  An/At.  After  this  date  the  R.  A.  E. 
values  were  at  irregular  intervals  and  the  Spacetrack  values  were  used.  As 
Fig.3  shows,  there  were  no  significant  differences  between  the  R. A.E.  and 
Spacetrack  values  when  both  were  available,  so  there  should  be  no  discontinu¬ 
ity  at  the  change-over  point. 


The  variation  of  perigee  height  for  a  polar  orbit  is  given  by 


p  U  Cn 

y  =  K  -  9*2  sin  u>  +  21  *4  sin  w  -  —  In  ~  km  , 


where  <x>  is  the  argument  of  perigee,  e  the  eccentricity,  and  H  (taken  as 

32  km)  the  density  scale  height  at  a  height  3H/2  above  perigee.  The  constant 

K  was  chosen  to  give  the  best  fit  to  observed  values  of  y  ;  the  values  used 

P 

were  176*0  km,  174*4  km  and  180*0  km  for  I966-5IA,  B  and  C  respectively. 

The  ’observed1  values  of  y  were  taken  from  R.  A.E.  orbits  for  51 A  and 

P  h 

B,  and  from  Spacetrack  bulletins  for  51 C.  The  3  R. A.E.  values  of  y^  and  the 

12  NASA  values  shown  in  Fig. 2  indicate  that  our  chosen  value  of  K  for  51 C  is 

probably  not  more  than  2  km  in  error.  It  should  be  noted  that  the  Spacetrack 

values  of  y^  in  Figs.1  and  2  require  a  correction  of  about  -7  sin  (0  km  to 

bring  them  into  agreement  with  the  usual  definition  of  perigee  height  (for 

further  explanation  see  Refs. 2  and  7).  In  Fig* 3  the  value  of  K  is  chosen  so 

that  the  curve  passes  through  the  group  of  values  between  MJD  39330  and  39330, 

5 

which  are  much  more  accurate  than  the  others  . 

2.3  Area/mass  ratios 

The  main  difference  between  the  three  orbits  arises  from  the  different 
area/mass  ratios  of  the  three  satellites,  giving  lifetimes  of  6,  13  and  9 
months  for  I966-5IA,  B  and  C  respectively.  The  orbital  acceleration  is  pro¬ 
portional  to  6  =  F  S Cp/m,  where  F  =  1  for  a  polar  orbit,  S  is  the  effec¬ 
tive  cross-sectional  area,  is  the  drag  coefficient,  and  m  is  the  mass 

of  the  satellite. 


For  1 966-51 B,  a  rectangular  box  measuring  33  x  28  x  23  cm,  the  effective 

2  2 

cross-sectional  area,  assuring  =  2*2,  was  between  0*078  m  and  0*115  m  , 


5 


2  2  2 

and  the  mean  value  of  S  was  taken  as  0*10  m  ,  giving  6  =  0*0129  m  /kg 
with  an  error  cf  ±1Q$  (s.d. ). 

The  shape  of  1 966-51 C  (erroneously  described  in  Ref. 4  as  'disc-like') 
is  that  of  a  regular  octahedron  of  side  0*23  m,  and  choosing  a  value  for  S 
posed  an  interesting  geometrical  problem.  The  four  principal  projections  of 
a  regular  octahedron  of  unit  side  are  illustrated  in  Figs.4a-d:  heavy  lines 
represent  edges  which  are  perpendicular  to  the  line  of  sight  and  therefore 
of  unit  length;  other  visible  edges  are  represented  by  continuous  light  lines, 
and  invisible  edges  by  broken  lines.  The  maximum  and  minimum  values  of  effec¬ 
tive  cross-sectional  area  are  shown  by  Figs.4&  and  4b,  giving  possible  extreme 
values  of  S  for  1 966-51 C  of  0*0529  m2  and  0*0374  m2  respectively,  if  the 
satellite  were  stabilized  in  either  of  these  directions.  It  is  more  likely 

that  the  satellite  was  spinning,  about  an  unknown  axis,  so  that  the  best  value 

2 

of  S  is  probably  that  given  by  Fig«4c,  i.e.  0*0458  m  ,  which  is  the  same  as 

for  spin  with  random  orientation,  as  well  as  being  close  to  the  mean  of  the 

2 

two  extremes  quoted  above.  We  therefore  took  S  =  0*0458  m  ,  with  an  error  of 
±8$  (s.d.).  With  =  2*2  and  m  =  5*0  kg,  this  gives  6  =  0*0202  m2/kg. 

As  the  size  and  weight  of  the  rocket,  1966-51  A,  were  not  accurately 

known  (it  exhibited  much  sharper  variations  in  brightness  than  a  normal  Agena 

rocket),  it  seemed  best  to  choose  a  value  of  6  to  give  values  of  density 

consistent  with  those  from  the  other  two  satellites.  The  value  of  6  used 
2 

was  0*031  m  /kg,  with  a  standard  deviation  which  will  be  assessed  in  section  3*2. 
2.4  Calculation  of  density 

Values  of  density  were  calculated  using  the  equations 


Pa  = 


0-157  (lOOOn)  /_e_Y  0-e)1/2 

.Jo  2  ,  l  aH*  /  T, 3/2 

10  n  6  \  /  (1+e)" 


1 


-  8e 


J2L\ 

1 6ae  / 


0*00335  cos  2<‘i  j 


•  •  • 


(2) 


(3) 


where  n  is  the  mean  angular  motion  in  rev/day,  a  is  the  semi  major  axis  in 
km,  and  H*  is  the  value  of  H  at  perigee;  y^  is  a  chosen  'mean'  perigee 
height  and  H'  is  the  value  of  H  at  a  height  y^  +  Equation  (2)  gives 

p  with  an  error  of  1$  (s.d.)  if,  as  here,  e  >  3H*/a.  For  all  three 

A 
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satellites  y^  was  chosen  as  177*5  ton*  while  H*  and  H'  were  taken  as  f 

27  km  and  30  km  respectively.  The  resulting  values  of  density*  pg  apply  to 
a  height  y^  +  =  191  km*  and  will  be  written  as  p^. 

3  RESULTS 

3*1  Values  of  density 

The  resulting  values  of  density  at  a  height  of  1 91  km,  ^ »  are 

shown  in  Fig. 5,  a  complicated  diagram  which  calls  for  extensive  explanation. 

Each  of  the  260  values  of  density  obtained  is  plotted  either  at  the 
epoch  of  the  orbit  determination  (for  the  R. A. E.  orbit),  or  at  the  mid-point 
of  the  time  interval  over  which  it  applies  (for  Spacetrack  values):  squares, 
circles  and  triangles  indicate  51  A,  B  and  C  respectively.  The  points  have 
been  joined,  by  unbroken  lines,  in  histogram  form:  since  the  Spacetrack 
values  of  n  =  An/ At  represent  an  average  over  the  time  interval  At,  the 
values  of  p ^  ^  are  also  averages  over  the  time  interval  At,  and  a  histo¬ 
gram  is  a  more  appropriate  presentation  than  joining  the  points.  The  values 
of  n  from  the  R.  A.E.  orbit  are  nominally  'instantaneous'  values,  at  the 
epoch  of  the  orbit  determination;  but  in  fact  the  variation  of  n  with  time 
was  taken  as  linear  in  most  of  the  orbit  determinations  of  1966-51  A,  so  that 
the  'instantaneous'  value  is  in  effect  a  mean  value,  over  a  time  interval 
whose  end  points  may  conveniently  be  taken  midway  between  the  successive 
epochs  of  the  orbit  determinations.  » 

Comparison  between  the  values  of  density  given  by  I960 -51  A,  B  and  C  in 
Fig. 5  is  difficult,  because  the  times  at  which  the  values  are  plotted  are 
different,  and  so  are  the  time  intervals  over  which  they  are  averaged.  In 
view  of  these  differences  and  the  'spiky'  nature  of  the  variations  in  density 
linked  with  geomagnetic  activity  -  the  peaks  at  MJD  39474  offer  a  clear 
example  -  comparisons  are  bound  to  be  inexact.  On  the  other  hand,  we  are 
unlikely  ever  to  have  any  better  sets  of  data,  unless  satellites  are  specially 
launched  for  this  purpose,  so  the  comparisons  are  worth  making. 

3.2  Cross-sectional  area  of  1966-51 A 

We  concentrate  first  on  comparing  the  mean  of  the  values  from  51 B  and  C 
with  the  values  from  51  A,  to  expose  any  large  variations  in  the  cross-sectional 
area  of  51  A.  Taking  the  mean  of  B  and  C  is  easier  said  than  done:  we  must 
average  the  values  of  p^^  over  time  intervals  that  are,  as  nearly  as  possible, 
the  same  for  both  satellites.  The  values  have  therefore  been  combined  into  , 
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groups,  as  shown  by  the  broken  lines  in  Fig. 5.  The  grouping  was  made  with  the 
aim  of  ensuring  that  the  end  points  of  the  corresponding  time  intervals  for 
51 B  and  C  should  if  possible  not  differ  by  more  than  12  hours,  and  should  never 
differ  by  more  than  24  hours.  This  grouping  of  the  values  obscures  some  of  the 
fine  detail:  but  such  a  loss  is  inevitable  if  a  valid  comparison  is  to  be  made; 
and  the  fine  detail  has  already  been  analysed  in  Ref. 2,  where  it  is  noted,  for 
example,  that  significant  increases  in  density  are  evident  at  the  times  of  17 
of  the  19  geomagnetic  disturbances  during  the  life  of  1 966-51 B. 

The  average  values  of  density  given  by  1 966-51 B  and  C  over  each  of  the 
chosen  time  intervals  are  written  between  the  two  graphs  in  Fig.5«  These 
values  sometimes  differ  by  as  much  as  3C$,  for  reasons  that  will  be  discussed 
later.  For  the  moment  we  assume  that  the  arithmetic  mean  of  the  values  from 
51 B  and  C  gives  the  best  indication  of  the  values  of  density,  and  these  mean 
values  are  superposed  as  a  dot-dash  histogram  on  the  values  of  density  from 
1966-51 A  at  the  top  of  Fig. 5.  Although  the  time  intervals  have  different  end¬ 
points  the  dot-dash  histogram  agrees  remarkably  well  with  the  values  of  p 
from  1966-51  A,  shown  by  the  histogram  with  unbroken  lines.  On  only  one 
occasion,  between  MJD  39400  and  39403,  does  the  average  of  the  values  from 
1 966-51 B  and  C  differ  by  more  than  1C$  from  the  average  value  given  by  1966-51A 
over  the  same  time  interval. 

This  unexpectedly  good  agreement  between  satellites  moving  quite 

independently  provides  a  conclusive  empirical  demonstration  that  (a)  the 

errors  inherent  in  the  method  are  less  than  1($,  (b)  the  effective  cross- 

sectional  area  of  1 966— 51 A  did  not  change  from  week  to  week  by  as  much  as 

2 

+1Q&,  and  (c)  the  assumed  value  of  6  for  51A,  0*031  m  /kg,  corresponding 
2 

to  S/m  =  0*0141  m^/kg,  was  correct  to  within  ±1($.  If  51A  had  been  a  normal 
Agena  rocket,  of  length  8  m  and  diameter  1*5  in,  the  expected  maximum  varia¬ 
tions  from  the  mean  cross-section  would  have  been  ±13$,  and  since  the  actual 
variations  are  smaller,  it  seems  very  probable  that  any  extra  equipment 
aboard,  responsible  for  the  satellite’s  rather  flashy  appearance,  did  not 
greatly  affect  the  cross-sectional  area;  i.e.  the  f3at  surfaces  causing  the 
flashes  were  probably  integral  with  the  body  rather  than  on  protruding  vanes. 

If  the  sise  and  shape  of  51 A  were  similar  to  a  normal  Agena,  the  mean  cross- 
sectional  area  would  be  about  10  m  and  the  mass  about  10/0*0141  ^  700  kg, 
which  is  close  to  the  usual  mass  of  an  empty  Agena  rocket.  Thus  the  flashing 
of  1966-51 A  was  probably  caused  by  fairly  trivial  alterations  to  the  basic 
design  of  the  Agena  rocket. 
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3. 3  Comparison  of  densities  from  1 966-51 B  and  51 C 

A  comparison  of  the  pairs  of  values  of  density  obtained  from  1 966-51 B 
and  C,  written  on  Fig.5>  shows  that  the  values  from  51 C  are  on  average 
greater  than  those  from  51 B.  The  possible  sources  of  this  bias  are  errors  in 
(a)  perigee  height,  (b)  cross-sectional  area  and  (c)  mass.  The  perigee  height 
of  1 966-51 B  should  be  accurate  to  rl  km;  the  perigee  height  of  51C  is  based 
on  rather  more  scattei’ed  observational  values  (see  Fig.  2),  and  the  standard 
deviation  is  estimated  as  1^  km,  corresponding  to  an  error  (s.d. )  of  3%  in 
p^.  The  estimated  errors  (s.d.)  in  the  cross-sectional  areas  of  1 966-51 B 
and  C  have  already  been  quoted  as  1C$  and  8j£  respectively.  Recent  data  on 
the  masses  of  the  satellites  suggest  that  the  values  used  should  not  be  in 
error  by  more  than  3$*  These  three  sources  of  error  are  enough  to  account 
for  the  16$  average  difference  between  the  densities  found  from  1966-51B  and 
C.  It  is  n  *ither  necessary  nor  possible  to  specify  how  much  each  of  the  three 
errors  contributes  to  the  total;  but  perhaps  the  most  likely  explanation  is 
that  the  area  taken  for  1966-51B  was  about  1C$  too  large,  and  that  the  remain¬ 
ing  ($>  difference  comes  from  errors  in  the  perigee  height,  cross-section  and 
mass  of  51C. 

3.4  Cross-sectional  areas  of  1 966-51 B  and  C 

The  variations  in  the  ratio  of  the  density  given  by  1  966-51 C,  p^^ 
say,  to  the  density  given  by  I966-5IB,  p^g  say,  provide  a  record  of  the 
variations  in  the  relative  cross-sectional  areas  of  the  satellites,  though 
we  have  to  remember  that  the  record  may  be  slightly  contaminated  by  error 
because  the  time  intervals  over  which  the  density  is  evaluated  are  not  exactly 
the  same  for  the  two  satellites.  We  write 


J*  = 


151L 


1  *1 6  p 


(4) 


51  B 


where  the  factor  1*1 6  has  been  inserted  so  that  p  has  a  mean  value  of  1; 

is  inversely  proportional  to  the  assumed  cross-sectional 


then,  since 


area 


351C’ 


P51C 

and  similarly  for  51B, 

C! 

|i  =  K 


513' 


J51C 


S51  (/“51  C 


£ 


f 


(5) 
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where  K  is  a  constant  which  it  is  convenient  to  write  in  terms  of  the  mean 
values  of  S^g  and  namely  S^g  and  respectively. 

The  41  values  of  p  obtained  from  the  values  of  ^  and  g  in 

Fig, 5  are  plotted  as  circles  in  Fig, 6.  The  variations  in  p  -  between  0*86 

and  1*14  -  are  smaller  than  might  be  expected.  For  if  Sc._  was  taken  1C % 

—  2 

too  large,  as  we  suspect,  vie  have  £L, _  =  0*090  m  ,  and  from  Ref. 2  the  limits 

_  2 

for  are  0*078  <  S^g  <  0*115  nr.  So,  if  =  0*0458  m  ,  as  assumed, 

and  0*0374  <  <  0*0 529  m^,  the  extreme  limits  for  p  are  0*75  <  p  <  1  *57* 

Since  in  fact  p  remains  betv/een  0*86  and  1  *14,  our  first  conclusion  from 
Fig. 6  is  that  it  provides  strong  observational  support  for  the  hypothesis  that 
the  cross-sectional  areas  of  spinning  satellites  very  rarely  reach  their 
extreme  values.  Indeed  Fig. 6  suggests  that  the  variations  in  cross-sectional 
area  are  much  less  than  the  estimated  bias  error  in  the  mean  area,  and  that 
the  s.d.  of  the  variations  may  be  taken  as  approximately  7  of  the  difference 
between  the  mean  and  extreme  values,  instead  of  about  ^ ,  as  we  usually  assume. 
Taking  the  factor  3-  would  give  variations  (s.d.)  of  1%  in  S^g  and  6%  in 
S^c,  i.e.  a  variation  of  9$  in  p:  the  actual  variation  (s.d.)  in  p  is 
8$,  so  the  factor  5  provides  an  adequate  margin  of  error. 

The  variations  in  p  in  Fig. 6  cannot  be  dismissed  as  mere  scatter; 
they  evidently  have  a  significant  structure.  To  remove  distracting  jaggedness, 
four  pairs  of  points  in  Fig. 6  have  been  averaged,  and  these  average  values, 

•  indicated  by  triangles,  have  been  used  when  joining  up  the  points  with  the 

unbroken  lines.  If  vie  accept  that  the  variations  in  p,  i.e.  in  the  ratio 
of  the  cross-sectional  areas  of  I966-5IB  and  C,  have  a  significant  structure, 
the  variation  in  p  will  be  regular  and  continuous  if,  as  is  likely,  the 
satellites  vie re  spinning  and  the  changes  in  cross-section  are  caused  by 
regular  precession  and  nutation.  A  tentative  curve  has  been  drawn  in  Fig. 6 
to  suggest  a  likely  form  for  the  variation.  Our  second  conclusion  from  Fig. 6 
is,  therefore,  that  it  reveals  a  regular  variation  in  S^g/S^g,  as 
indicated  approximately  by  the  broken  curve. 

The  third  conclusion  to  be  drawn  from  Fig. 6  is  concerned  with  the 
remarkably  small  scatter  of  the  points  about  the  broken  curve,  about  2fo  rms. 

It  must  be  remembered  that  the  values  plotted  in  Fig. 6  represent  the  ratios 
of  values  of  density  obtained  at  slightly  different  times,  averaged  over 
different  time  intervals,  from  satellites  moving  in  independent  orbits  and 
having  variable  cross-sectional  areas.  Since  the  corresponding  time  intervals 
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for  51 B  and  C  often  differ  by  as  much  as  2C$,  and  the  change  in  density  between 
one  interval  and  the  next  is  also  often  as  much  as  2C$,  errors  in  p  of 
could  arise  from  this  source  alone;  since  the  actual  scatter  in  p  is  only 
about  2?o,  errors  inherent  in  the  method  of  density  determination  must  be  less 
than  2?o. 

To  illustrate  the  contribution  of  differences  in  the  time  intervals  to 
the  scatter,  consider  the  point  furthest  from  the  broken  curve  in  Fig. 6,  that 
at  MJD  39470.  As  Fig. 5  shows,  the  value  of  p^c  at  this  point  is  slightly 
'contaminated'  by  the  beginning  of  the  magnetic  storm  at  MJD  39478,  whereas 
the  effect  of  this  storm  on  is  concentrated  in  the  subsequent  value  of 

I  ^  ft  - — 

density.  Averaging  of  the  two  successive  values  of  p  should  noticeably 
reduce  the  scatter  -  and  does  so,  as  shown  by  the  triangle  in  Fig. 6  at  day 
39472.  Examination  of  other  points  not  very  near  the  curve  leads  to  similar 
conclusions. 


3.5  Discussion  of  variations  in  u 


In  considering  the  variations  of  p,  the  latitude  of  perigee  is  relevant. 
The  scales  at  the  top  of  Fig.  6  show  the  values  of  u>,  the  argument  of  perigee 
for  51 B,  and  the  latitude  of  perigee,  <j>  »  The  variation  of  p  shows  a 
consistent  pattern:  the  four  maxima  in  the  broken  curve  occur  at  0  +70° 

and  the  three  minima  occur  when  perigee  is  near  the  equator  (<j>  cl  O). 


Details  of  the  spin  behaviour  of  the  two  satellites  are  not  known,  but 
it  is  believed  that  1966-51C  was  spinning  at  rates  close  to  its  design  value 
of  12  rev/min.  If  so,  its  cross-sectional  area  should  not  change  by  more  than 
a  few  per  cent,  and  the  variations  in  p  would  be  primarily  due  to  variations 

m  g. 

Certainly  1 966-51 B,  a  rectangular  box  measuring  33  x  28  x  23  cm,  had 

more  opportunities  to  change  its  cross-sectional  area.  The  satellite  was 

designed  to  be  magnetically  stabilized;  but  it  entered  an  incorrect  orbit 

0 

on  which  the  aerodynamic  forces  at  perigee  were  about  10  times  greater  than 
on  the  intended  orbit.  On  the  actual  orbit  the  magnetic  stabilization  is 
unlikely  to  have  been  fully  successful;  but  it  may  have  been  partially 
effective,  for  the  limited  evidence  from  radar,  radio  and  optical  tracking  is 
all  consistent  with  a  slow  rotation  rate,  which  may  have  been  an  angular 
oscillation  rather  than  true  spin.  If  so,  the  minimum  cross-sectional  area 
might  have  occurred  when  perigee  was  near  the  equator,  when  the  longest  side 
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(33  cm)  would  be  horizontal  if  it  was  aligned  with  the  magnetic  field;  SR. _ 

2  2  ■>'a 
could  then  be  as  low  as  0*078  m  ,  as  compared  with  the  probable  average  value 

2 

(see  section  3*4)  of  0*090  m  .  This  may  explain  the  minima  in  |i  when  perigee 
is  near  the  equator.  The  maxima  in  S^g  at  <f>  ^  i 70 0  cannot  so  easily  be 
accounted  for:  the  north  and  south  magnetic  poles  are  at  latitudes  of  74°N 
and  68°S;  if  they  are  responsible  for  the  maxima,  however,  there  should  also 
be  maxima  at  w  =  110°  and  250°,  when  perigee  latitude  is  70°.  It  is  true 
that  there  are  maxima  in  (i  at  MJD  39317  (w  -  100°)  and  MJD  39401  (to  =  260°), 
though  they  might  not  be  recognized  except  by  those  who  were  looking  for  them. 
Despite  its  deficiencies,  this  'magnetic  explanation'  may  well  be  correct. 

If  the  variations  in  (i  are  caused  by  variations  in  the  cross-sectional 
area  of  1966-51B,  the  values  of  air  density  obtained  from  that  satellite  should 
be  corrected  to  allow  for  the  variations  in  S^g.  If  this  correction  was  made, 
the  maxima  and  minima  of  the  semi-annual  variation  in  density  (averaged  over 
30  days)  would  be  as  follows: 

Af^c  T  -1--  3*7  x  10_1°  kg/m^ 


1 966  July  minimum 

1966  October  maximum 

1 967  January  minimum 
1 967  March  maximum 


(unchanged) 


In  fact  we  cannot  be  sure  that  this  full  correction  should  be  made,  and  the 
safest  procedure  is  to  take  the  mean  of  these  and  the  original  values:  the 
amplitude  of  the  semi-annual  •  jriation  in  density  would  then  be  reduced  from 
the  value  of  1  *45  given  in  Ref* 2  to  1 *40. 

We  also  suspect  (see  section  3*4)  that  the  mean  value  of  S^g  should 
be  rather  smaller  than  was  assumed  in  Ref. 2.  If  so,  all  the  values  of  density 
in  Ref. 2  should  be  slightly  increased.  We  cannot  specify  the  increase  exactly, 
but  it  is  unlikely  to  exceed  1C$. 

4  CONCLUSIONS 

(l)  The  effective  cross-sectional  area  of  the  Agena  rocket  I966-5IA 
changed  by  less  than  1C$  from  week  to  week,  and  the  average  area/mass  was 
0*014  m  /kg.  The  size  and  mass  were  probably  similar  to  a  normal  Agena,  and 
the  rather  flashy  appearance  was  probably  caused  by.  fairly  minor  alterations 
to  the  basic  shape. 
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(2)  The  quantity  p  plotted  in  Fig. 6  represents  the  (normalized)  ratio 
of  the  effective  cross-sectional  area  of  1 966-51 B  and  51 C,  i.  e.  S^g/S^, 

as  deduced  from  the  variations  in  orbital  decay  rates  for  the  two  satellites. 
Although  theoretically  capable  of  ranging  between  0*75  and  1 *57,  p  remains 
between  0*86  and  1 *14,  thus  confirming  the  usual  assumption  that  the  effective 
cross-sectional  areas  of  satellites  rarely  approach  their  extreme  values. 

(3)  Fig. 6  shows  that  p  undergoes  regular  oscillations  with  an 

amplitude  of  about  1C$.  This  should  probably  be  interpreted  as  a  variation 
in  since  S,.^  is  likely  to  have  remained  fairly  constant.  If  so, 

Sj-.jg  had  minimum  values  when  perigee  was  at  the  equator  and  maximum  values 
at  latitudes  near  70°.  The  most  likely  explanation  is  that  the  intended 
magnetic  stabilization  of  1 966-51 B,  a  rectangular  box,  was  partially  effective 
and  that  the  maxima  and  minima  in  S^g  occur  near  the  magnetic  equator  and 
poles. 

(4)  If  S^g  varied  in  the  manner  indicated  by  p,  some  small  amend¬ 
ments  should  be  made  to  the  values  of  air  density  deduced  in  Ref. 2  (in  which 
S^g  was  assumed  constant).  In  particular  the  amplitude  of  the  semi-annual 
variation  in  air  density  in  1 966-7  at  a  height  of  1 90  km  should  be  changed 
from  1  -45  (as  given  in  Ref. 2)  to  1 *40. 

(5)  The  values  of  p  in  Fig. 6  are  obtained  from  the  ratio  of  orbital 
decay  rates  of  1966-513  and  C,  evaluated  over  slightly  different  time  intervals, 
and  errors  of  tjfo  (s.d.)  would  be  expected  from  this  source,  in  addition  to  any 
errors  due  to  inaccuracies  in  orbital  data  and  the  theory  used.  In  fact  the 
scatter  of  the  values  of  p  about  the  curve  is  2$  (s.d.):  this  severe  test 
provides  conclusive  empirical  confirmation  that  the  errors  inherent  in  the  data 
and  method  of  analysis  give  rise  to  errors  of  less  than  2$. 
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Fig.  3  Orbital  decay  rate,  ri, and  perigee  height, 
y.p,  for  1966-  51 A  (Secor  6  rocket) 


Fig.  4  a-d 
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(A—  I ) 


(b)  Along  plan*  of  symmetry 
(X*l JJ2  -  0-707) 


(c)  Perpendicular  to  face 
(W 3/2  -0  -866) 


(d)  Parallel  to  face 
(A  -72/v/3  -0-816) 


Fig4a-d  Visible  areas  A  of  a  regular  octahedron  with  edges 
of  unit  length,  along  four  principal  lines  of  sight 


